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Abstract 

Wc show that the problem of counting the number of Eulerian circuits in an undirected 
graph is complete for the class #P. 

1 Introduction 

Every basic text in graph theory contains the story of Euler and the Konigsberg bridges, 
together with the theorem that guarantees the existence of a circuit traversing every edge of a 
graph exactly once, if and only if the graph is connected and all vertices have even degree. 

Later in the text one might find the "Matrix Tree Theorem", which provides an efficient 
algorithm for counting the number of spanning trees of a graph; and an application of this 
and the so-called "BEST" Theorem (see below) to count the number of Eulerian circuits in a 
directed graph. But what about counting Eulerian circuits in an undirected graph? 

This problem is clearly in the class #P (introduced by Valiant in the 1970s), since 
it is easy to check whether a candidate circuit traverses each edge once. Since there is no 
known efficient way to count Eulerian circuits, it is natural to suspect that the problem is #P- 
complete, and thus presumably very difficult — especially in view of Toda's result [§], which 
implies that one call to a $P oracle suffices to solve any problem in the polynomial hierarchy 
in deterministic polynomial time. 

Researchers have shown a myriad of graphical counting problems to be #P-complete, in- 
cluding Hamilton circuits ^2], acyclic orientations [§], and Eulerian orientations 0. Yet, the 

complexity of counting Eulerian circuits remained open (see e.g. [3j, Open Problem on p. 5, at 
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the end of Section 1.1) for 25 years — rather mysteriously, especially considering the simplicity 
of the reduction below. 

Our approach is to show that, with the help of an oracle which counts Eulerian circuits, a 
Turing machine can count the number of Eulerian orientations of any given graph in polynomial 
time. The latter problem was shown to be #P-hard by Mihail and Winkler [7J |S] and for 
completeness we provide a second proof below. 

Both reductions proceed by enumeration modulo various primes, a technique introduced 
originally by Valiant and utilized later by the authors [2] to settle another long-open 
complexity problem, counting the linear extensions of a partially ordered set. 

To show later that we can reconstruct a number uniquely from its values modulo a small 
set of primes, it is useful to have a technical lemma such as the following. 

Lemma 1.1. For any n > 4, the product of the set of primes strictly between n and n 2 is at 
least n\2 n . 

Proof. We use some facts from Hardy and Wright [I], Chapter 22, concerning the functions 

log n I log 2 

i?(n) = logn p <n^' wnere P runs over all primes less than n, and ijj{n) = r d{n 1 ^ % ). 

i=l 

From |3] we find that i9(n) < 2nlog2 for n > 1, and that ip(n) > |nlog2 for n > 2. 
We are interested in the quantity V = i9(n 2 ) — i?(n). From the above facts, we have: 

> in 2 log2- •2nlog2-2nlog2 
> nlogn > log(n!2 n ) , 

at least provided n > 150. The inequality for 4 < n < 150 is easily verified by direct calculation. 

□ 

It is evident that this lemma is not tight: it is possible to replace the n 2 upper limit by 
iTnlogn, for some suitably large K. 

2 Circuits, Orientations, Arborescences and Orbs 

For us a graph G = {V, E) will be finite and undirected, with no loops or multiple edges; if 
multiple edges are permitted we use the term multigraph. A circuit C of G is a closed path, 
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with a direction but no distinguished starting point; it is Eulerian if it traverses each e € E 
exactly once. Of course the possession of even one Eulerian circuit implies all degrees are even, 
and it will be convenient for us to denote the degree of a vertex v by 2d v instead of d v . 

An Eulerian orientation of G is an orientation of its edges with the property that each 
vertex has the same number (namely, d v ) of incoming and outgoing arcs. Any Eulerian circuit 
induces an Eulerian orientation by orienting each edge in accordance with its direction of 
traversal. 

If a particular starting edge is chosen for the Eulerian circuit C, originating say at vertex 
r, then C also induces a spanning tree T = {exit(f) : v ^ r} where exit(v) is the last edge 
incident to v used by C before its final return to r. When oriented according to C, T becomes 
an in-bound spanning tree, or arborescence, rooted at r. 

Let us fix a root r G V and denote by the term orb a pair consisting of an Eulerian 
orientation and an arborescence (for that orientation) rooted at r. If an orb is specified, it 
is easy to construct a corresponding Eulerian circuit C: simply begin walking from the root 
r, following any unused outgoing arc from each vertex v, except that the tree arc exiting v 
is avoided as long as possible. Since there are (d v — 1)1 ways to order the non-tree outgoing 
arcs from v, and d r \ from r, the number of ways to construct C is precisely d r \ Y\ v ^ r {d v — 1)!. 
However, this over-counts Eulerian circuits (as we have defined them) by a factor of d r since each 
circuit passes d r times through r. Hence orbs and Eulerian circuits are in perfect n^ey^ — -0'" 
to-1 correspondence, and thus counting Eulerian circuits is equivalent to counting orbs. This 
result is sometimes known as the "BEST" Theorem after de Bruijn, van Aardenne-Ehrenfest, 
Smith and Tutte, although the former two should perhaps get additional credit as the original 
discoverers. 

Note that for any particular orientation, one can use the BEST Theorem together with the 
Matrix- Tree Theorem jlUj (due to Tutte, but implicit in Kirchhoff 's work of 1847) to count 
Eulerian circuits. The difficulty in using this approach to count (or even approximate) the 
number of Eulerian circuits in G is that the number of arborescences in an orientation can 
vary enormously. 

Theorem 2.1. Counting Eulerian circuits is #P- complete. 

Proof. It suffices to reduce the problem of counting Eulerian orientations to counting orbs in a 
multigraph; the latter is equivalent to counting orbs in a simple graph since multiple edges can 
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be subdivided without affecting the number of orbs. Let us therefore assume that an Eulerian 
graph G = (V = {1, . . . , n}, E) has been given, and that we have an oracle for counting orbs 
in any multigraph. We wish to compute the number N of Eulerian orientations of G. 

We construct for any odd prime p, a graph G p whose number of orbs is equivalent to ./V 
modulo p. The construction is shockingly simple: each edge e of E is replaced by p parallel 
edges ei, . . . , e p , and a new node is added which is adjacent to every node v of V by two 
parallel edges, e(v) and e'(v). We take as the root of all orbs. 

Figure Q shows a small G and the resulting G p when p = 3. 




Figure 1: Reducing eulerian orientation count to ORB COUNT 

In G p , the type r of an orb is a function from E to {0, 1, . . . ,p} x {T, F} which tells how 
many of e±, . . . , e p are oriented from the smaller to the larger-numbered vertex, and whether 
any is a tree edge ("T" means "yes", "F" means "no"). A type r is special if r(e) £ {0,p} x {F} 
for every e££. 

If an orb belongs to a special type, then the common direction of the parallel edges corre- 
sponding to each eGfi provides an orientation of G. This orientation is Eulerian, as otherwise 
the in-degree and out-degree of some vertex of G p will differ by at least 2p — 2. Therefore, in 
an orb of special type, for each vertex v of G p , exactly one of e(v) and e'(v) is directed towards 
the root 0, and this is the arc that carries the edge of the associated tree directed away from 
v. Thus, the number of orbs of special type is precisely 2 n x N. 

On the other hand, we claim that the number of orbs of any non-special type r is a multiple 
of p. To see this let e be such that r(e) = (k, X) ^ {0,p} x {F}. Suppose first that < k < p 
and X = F; then the orbs of type r can be partitioned into (?) equal parts according to which 
of oriented from the smaller to the larger-numbered vertex, and of course (|) is 
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a multiple of p. 

If < k < p and X = T, then the part sizes are multiples either of k(^j or (p— k)(f) since 
we must also decide which of the correctly-oriented arcs belongs to the tree. 

Finally, if k £ {0,p}, and X = T, we partition according to which of the now-parallel arcs 
is in the tree and there are p choices. 

It follows that the total number of orbs is equivalent to 2 n x N modulo p, and thus we 
can compute N mod p. We repeat this process for every prime p between m and 2m, where 
\E\ = m, and apply the Chinese Remainder Theorem to nail N. □ 

A similar argument, equally simple, can be used to reduce not-all-equal 3-SAT COUNT 
(shown in [Hj to be #P-complete) to counting eulerian orientations. Given an instance 
of not-all-equal 3-SAT count, the graph is provided with a vertex for each literal and 
another for each clause, plus one spare vertex s. Each clause is connected by a single edge to 
its three literals and s. Each literal vertex is given p parallel edges to its mate (where p is 
a prime larger than the number of appearances of any literal), and a total of p other edges, 
consisting of some number (as mentioned above) to the clauses in which it appears, with the 
remainder going to s. 

Figure El shows the construction for a particular 2-clause, 3- variable instance with p = 3. 



X L/ Z X Lj Z 




x x y L) z z 

Figure 2: Reducing not-all-equal 3sat count to eulerian orientation count 

As in the proof above, orientations which fail to align in parallel all p edges associated with 
any given variable fall into classes of size mod p. Each of the other "special" orientations has 
the property that for each variable x, either every edge from x to a clause vertex points outward 
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and every edge from x to a clause vertex points inward, or vice- versa. Such an orientation 
corresponds to a satisfying assignment, since we cannot have three of the four edges incident 
to a clause vertex pointing the same way. Conversely, given a satisfying assignment, we orient 
all the edges between x and x towards the true literal, orient all other edges incident with a 
literal vertex away from true literals and towards false literals, and orient the edge between 
each clause vertex and s in the necessary manner. This ensures that the in-degree is equal 
to the out-degree at every vertex other than s, and therefore also at s. Hence there is an 
exact correspondence between special orientations and satisfying assignments, and we proceed 
as before. 
Remark 

Still open is the question of whether there is a fully polynomial randomized approximation 
scheme ("fpras") for counting Eulerian orientations (as there is, e.g., for Eulerian orientations 

[H] ). We believe that there is, and even that a particular Markov chain whose states are orbs 
and near-orbs mixes rapidly. We hope and expect that this question will not remain open for 
another 25 years. 
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